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Abstract. The s2s-OVCA is a cellular automaton (CA) hybrid of the optimal velocity 
(OV) model and the slow-to-start (s2s) model, which is introduced in the framework of the 
ultradiscretization method. Inverse ultradiscretization as well as the time continuous limit, 
which lead the s2s-OVCA to an integral-differential equation, are presented. Several traffic 
phases such as a free flow as well as slow flows corresponding to multiple metastable states 
are observed in the flow-density relations of the s2s-OVCA. Based on the properties of the 
stationary flow of the s2s-OVCA, the formulas for the flow-density relations are derived. 
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1 Introduction 

Self-driven many-particle systems have provided a good microscopic point of view on the vehicle 
traffic [3, 5]. The optimal velocity model [1] gives a description of such a system with a set 
of ordinary differential equations (ODE). It is a car-following model describing an adaptation 
to the optimal velocity that depends on the distance from the vehicle ahead. Another way of 
describing such systems is provided by cellular automata (CA). For example, the elementary 
CA of Rule 184 (ECA184) [16], the Fukui-Ishibashi (FI) model [4] and the slow-to-start (s2s) 
model [12] are CA describing vehicle traffic as self-driven many-particle systems. 

Studies of the self-driven many-particle systems have been wanting a framework that com¬ 
mands a bird’s eye view of both ODE and CA models in a unified manner. Ultradiscretiza¬ 
tion [14], which gives a link between the Korteweg-de Vries (KdV) equation and integrable soli- 
ton CA [11], is expected to provide such a framework, for it can be applied to non-integrable sys¬ 
tems, too. As a first step toward such a framework, an ultradiscretization of the OV model [10] 
was presented. A specific choice of the OV function enabled the ultradiscretization of the OV 
model without any other specialization. We should note that another ultradiscretization of the 
OV model with essentially the same OV function as above was derived from the modified KdV 
(mKdV) equation, which is an effective theory around the critical point, with specializing its 
solutions to traveling wave solutions [7] . The latter ultradiscretization of the OV model depends 
on the ultradiscretization of the mKdV equation, which is an integrable soliton equation with 
rich accumulation of the studies of integrable discretization and ultradiscretization. The former 

*This paper is a contribution to the Special Issue on Exact Solvability and Symmetry Avatars in honour of 
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one, on the other hand, has nothing to do with integrability, which indicates a possibility to 
expand the scope of the ultradiscretization beyond the integrable models. Thus the two ultra¬ 
discretizations make a clear contrast regarding to integrability. An early search for a CA-type 
OV model dates back to 1999 [5, 6], which was done from a phenomenological point of view to 
highway traffic. 

The former ultradiscretization of the OV model [10] provided a foundation toward a hybri¬ 
dization of the OV model and the s2s model, and it lead to the s2s-OVCA [9] indeed, without 
the aid of the integrable models. The s2s-OVCA is a CA-type hybrid of the OV model and 
the s2s model. As we shall see, the equation of the s2s-OVCA generally involves three or more 
times, or higher order time-differences, in other words. As far as the authors know, difference 
equations involving higher order time-differences yet want a thorough study from a point of 
view of the ultradiscretization. Besides an interest from the traffic theoretical point of view, an 
interest toward a new horizon of the scope of the ultradiscretization motivates us to introduce 
and study the s2s-OVCA. As we shall show in Section 2, the s2s-OVCA reduces to an ODE that 
is an extension of the OV model in the inverse-ultradiscrete and the time-continuous limits. 

It was observed by numerical experiments that motion of the vehicles described by the s2s- 
OVCA went to stationary flow in the long run, irrespectively of the initial configuration [9, 15]. 
It was also observed by numerical experiments that the flow-density relation for the stationary 
flow of the s2s-OVCA was piecewise linear and flipped-A shaped with several metastable slow 
branches [9] . Exact expression for the flow-density relation was given by a set of exact solutions 
giving stationary flows of the s2s-OVCA [15]. The flipped-A shaped diagram captures the 
characteristic of observed flow-density relations [3, 5]. Some other CA type models that gave 
a flipped-A shaped flow-density relation with a metastable branch was also reported [2], We 
shall explain in Section 3 the flow-density relation of the s2s-OVCA based on the properties of 
the stationary flow which was numerically observed [9]. 


2 s2s-OVCA and its inverse ultradiscretization 


The s2s-OVCA is given by a set of difference equations below 


_L1 /no 

x k =4 + min(mm(x;f + " 



( 1 ) 


where the integers no > 0, vq > 0 and x^, k = 1,2,..., K, are the monitoring period, the top 

speed and the position of the car k at the n-th discrete time. Note that the definition of the 

N 

symbol min is 
k =o 


N 

min(afc) := min(ao, aq, < 22 ,..., ajv). 
k =0 

The equation (1) is called an ultra-discrete equation in the sense that it is a difference equation 
which is piecewise linear with respect to the dependent variables x JJ. The s2s-OVCA includes 
the ECA184 (no = 0, uo = 1) [16], the FI model (no = 0) [4] and the s2s model (no = 1, 
vq = 1) [12] as its special cases. 

Since the second term in the right hand side of equation (1) gives the speed of the car k 
at the time n, the s2s-OVCA describes many cars running on a single lane highway in one 
direction, which is driven by cautious drivers requiring enough headway to go on at least for no 
time steps before they accelerate their cars. The equation (1) also means that the car slows 
down immediately when its headway becomes less than its velocity. Thus the monitoring pe¬ 
riod no describes asymmetry between acceleration and deceleration of the cars. It is said that 
the acceleration times are about five to ten times larger than the braking times [5]. 
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Without loss of generality, we can assume that the cars are arrayed in numerical order, 
x5 < x® < ■ ■ ■ < x® K , which is also assumed throughout below. Then the number of empty cells 
between the cars k and k + 1 for any k is always non-negative, i.e., 


Ax?-1>0. 


( 2 ) 


It is obvious that the inequality holds for n = 0. We assume that the inequality holds up to 
some n, as the induction hypothesis. The induction hypothesis as well as the definition of min 
assure the inequality 


( n 0 / \ 

min(Ax? -n — l),xoj < Ax? — 1 
for any k. Using equation (1), we get an expression of Ax? as 

Ax? + — 1 = Ax? — 1 + min^ruin(Ax? - ? — l), xo j — min?min(Ax? - ” — l), voj 

( n 0 , \ 

min(Ax?“? - l),x 0 ) + 


( 3 ) 


/ n 9 

Ax? — 1 — mini min 

\n '=o 


■ 0) 


The inequality (3) and the equation (4) show that the inequality (2) holds for n + 1. The 
inequality (2) means that both overtake and clash are prohibited by the s2s-OVCA. 

We should note that the s2s-OVCA is obtained from a difference equation by a limiting proce¬ 
dure named ultradiscretization [14], which generates a piecewise-linear equation from a difference 
equation via the limit formula 


/ N \ 

lim <5xlog( bke ak ^ Sx ) = max(ao, ai, a, 2 , ■ ■ ■, oat) =: max(ofc), 

c-H-0 \ z —' / ' ' k= 0 

v fc =0 7 


Sx —>-+0 

where arbitrary numbers bk must be positive. The equation (5) is rewritten as 

N ' -i JV 

= min(a fc ), 

k =0 


( 5 ) 


lim 5xlog{y"b k e- ak / Sx 
Sx —>-+0 ' ^ 


k =0 


for min(ao, a±, a 2 , ■ ■ ■, ajv) = — max(—ao, —ai, — < 22 ,..., — citv). 

For the sake of convenience in the calculation below, we introduce two parameters, xo and 5t, 
in the s2s-OVCA 


" n+1 = xl + min (mm (Ax? ”') - x 0 ,v 0 St) =: x? + x? pt (A eff x?)(5t, 


\n'=0 


( 6 ) 


where A e ffX? := min?? =0 (Ax?“” ). The parameters xq and St are the length of a cell, which 
corresponds to the space occupied by a single car or the length of a car itself in the shortest 
case imaginable, and the discrete time-step, respectively. For simulation of highway traffics, 
the length of a cell xo and the discrete time-step St are usually chosen as xo = 7.5 m and 
xo = 5 x || [5]. But we do not consider specific values of the parameters in the s2s-OVCA and 
regard them as generic. The two parameters xo and St were set to be unity in equation (1). 
Introduction of xq into the inequality (3) gives 


Ax? — xo > 0 


( 7 ) 


for any n and k , which is shown by induction with the aid of an inequality 

( no / \ 

min(Ax? -n — xo),xo(5fJ < Ax? — xq 
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and an expression of Ax^J derived from equation (6) 

/ n o / \ r / / 

Ax^ + — xq = min( min(Ax^” — xo),^oJ + AxJ? — xq — minf min(Ax^ _n — xq),uo 


that correspond equations (3) and (4), respectively. 

Since we have the inequality (7) for the headway Ax£ — xo, the effective headway A e ffx£ — xo 
is also always non-negative, A e ffxj? — xo > 0, for any k. With the aid of the identity 

min (A B) = A — max(0, A — B) = max(0, A) — max(0, A — B) 

for any A > 0, the optimal velocity function v^ pt (x)5t := min(x — xo,xo 5t) in the s2s-OVCA is 
expressed as 

t>o pt (x)<5t = max(0, x — xo) — max(0, x — xo — vq St), 
for any x > 0. It is given by the ultradiscrete limit 5x —> +0 of a function 


Vo P t(x)8t = dxlog 


l-|- e —=^0 / <5a: 

-^_^_Q(x — x Q—VQSt)/Sx 

l+ e -(a:o+' u O <5t )/ ,5:r 


which is an inverse-ultradiscretization of the optimal velocity function ^opf Note that we have 
introduced arbitrary coefficients so as to make ^opt (0) = 0. In a similar way to the above 
calculation, an inverse-ultradiscretization of the effective interval A“ ff x^ is also obtained as 


A efAfc : = <^log 


( n 0 

y -—-— 

is n » + 1 


-l 


Therefore an inverse-ultradiscretization of the us2s-OVCA is given by x£ +1 = x^+v^ pt (A ^ {{ x^)5t, 
which is explicitly written as 


xl +l = x'l + 5x1 log 


, ^ e -(A xl~ n -x 0 )/Sx 

1 + (E- 

n '=0 


n 0 + 


- log 


>+ E 


_ n °, e -(Ai^ n -x 0 -v 0 St)/Sx \-! 


n '=0 


n 0 + 1 


-xo)/dx \ —in 

- -J - log(l + e-*°/ fa ) 

+ log(l + e- (a:o+,Jo ' 5t) / fa )|. 


( 8 ) 


In other words, the s2s-OVCA is given by the ultradiscrete limit 5x —>• +0 of the above difference 
equation (8). 

Since equation (8) is rewritten as 


n+1 


5t 


-H-s" 0 * 


n ° e -(A x n k - n -XO -v 0 St)/Sx 

1 + (E- 

■n '=0 


n 0 + 1 


-i- 


- log 


. e -(A X n k ~ n -xo)/Sx 

1 + ( E- 

n '=0 


— 1- 


+ 


= Vo 


n 0 + 1 

log (1 + e -(*o+v 0 5t)/6x ) + + e -x 0 /5x ) 


i+E 


5t 


22. e -( A -r n 


n '—0 


n 0 + 1 


) -v 0 {l+ e X0 / Sx y 1 + O(5t), 
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the above difference equation (8) goes to an integral-differential equation in the continuum limit 
5t —> 0 as follows 

e —(A x h {t-t')-x 0 )/Sx^£ 

A ^^+1 _ ryfl 

where to : = no5t, -#■ = lim k x . —- and 
’ dt St—>o 5t 


d Xk 
dt 


= v Q 


l+ hk 


to 


) -v 0 {l+e x °/ 5x ) 


-l 


(9) 


"0 —(Ax?~ n -x 0 )/5x 1 /-to 

lim V - -= - [ e -(^k(t-t')-x 0 )/Sx dt i 

n '=0 


St—^0 ^ TT-o ~"t"~ 1 to JO 


In terms of an optimal velocity function and an effective distance 
v op t{x) :=v 0 ‘ 


AeffX' fc (t) := Jxlog 


\ _)_ q-(x—xo)/Sx y -)- qXq/Sx 

to \ — 1 

e -Ax fc (t-t')/5x^/ 


'0 


the above integral-differential equation is expressed as 
= ^opt (AgffaJfc(t)) • 

Since the effective distance A e ffXfc(t) goes to Ax k [t) in the limit below 

1 


A Xk{t — to) = lim log 
h—yto 


r ° e -A x k {t-t’)/5x d £ 


Jo ~ h J h 

this integral-differential equation is an extension of the Newell model [8] 
= v op t(Ax k (t - t 0 )), 


( 10 ) 


which is a car-following model dealing with retarded adaptation to the optimal velocity deter¬ 
mined by the headway in the past. 

Replacement of t with t+to in equation (10) and the Taylor expansion of Xfc(f+to) = Vk{t+to) 
yield 


, . ... . . . . dv k 1 d?Vh 9 

Uopt(Ax fc (f)) = v k (t + to) = Ufc(t) + ■ to + ^ ~ l -’ 

which is equivalent to 

lit + YW' tG + '" = ^K P t(AXfc(t)) -Ufc(t)). (11) 

The equation of motion of the OV model 

-jf = ^(«qpt(AXfc(t)) - Ufc(t)) 

is given by neglecting the higher order terms in the left hand side of the equation (11). 

The discussion shown above in this section shows how the inverse ultradiscretization and the 
time continuous limit connect the s2s model and the Newell model, which approximates the OV 
model, through the s2s-OVCA. 










6 


H. Ujino and T. Yajima 


J3 

<u 

> 


a 

o 


c/l 

o 

Oh 



Time 



Figure 1. 


The spatio-temporal pattern (left) and the flow-density relation (right) of the s2s-OVCA [9]. 


3 Flow-density relation 

Fig. 1 gives typical examples of the spatio-temporal pattern showing jams and the flow-density 
relation of the s2s-OVCA [9]. In the numerical calculation, the periodic boundary condition is 
imposed and the length of the circuit L , which is the same as the number of all the cells, is fixed 
at L = 100. The maximum velocity vq and the monitoring period no are vq = 3 and no = 2. 
The number of the cars K in the spatio-temporal pattern is set at K = 30. 

The spatio temporal pattern shows the trajectories of the cars. As we can see, irregular 
motion of cars is observed in the early stage of the time evolution, 0 < n < 30, where n is the 
time. But after that, the flow of the cars become stationary in the sense that length of the jam 
is almost constant and that cars with intermediate speeds appear only temporarily. 

The flows Q in the flow-density relation are computed by averaging over the time period 
800 = ni < n < ?rf = 1000, 


Q: = 


K n f 






n .n ._ n +1 _ n 
u k ■— - b k x ki 


( 12 ) 


in which the traffic is expected to be stationary in the above mentioned sense. The car density p 
is given by p := As we have mentioned before, the flow-density relation of the s2s-OVCA is 
piecewise linear and flipped-A shaped with several metastable slow branches. 

The flow-density relation shown above is derived by admitting the features of the flow of 
the s2s-OVCA. Namely, the flow of the s2s-OVCA goes to one of the stationary flows in the 
long run. The stationary flows consist of the free flow in which all the cars run at the top 
speed vo and the slow flows that always contain slow cars running at the minimum speed u“ in , 
0 < u“ n < Vo, which remains constant. Formation of the line of slow cars corresponds to that of 
traffic jam. In the slow flows, lengths of the jams are almost constant and fluctuate periodically. 
Our previous paper [15] gives a set of such stationary flows. 

First we shall deal with the free flow and its flow-density relation. Since all the cars run at the 
top speed, vj? = vo, V k , all the headways respectively remain constant, Ax£ +1 = Ax ^ > vo + 1, 
Vk. Hence the flow Q is also constant in the future. Using the definition of the flow (12) with 
n\ = rif = n, we have 


Q 



k =1 


PV 0 , 


which gives the straight line in the flow-density relation with a positive inclination that is equal 
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Figure 2. The formula for the flow-density relation of the s2s-OVCA [15]. 

to the top speed vq. For example in Fig. 2, the flow-density relation of the free flow which is 
labeled with v = 3 is Q = 3 p, since vo = 3 in this case. 

Next we shall consider the slow flows and their flow-density relations. Let us see a specific 
solution of the s2s-OVCA starting from the following initial configuration 

0: 1u2u3uuu4uuuuuuu5ulj6u7u8lj9uuuuuuu0ljuuu • (13) 

Note that the number 0 at the leftmost shows the time. The digits and the blank symbols u 
in the above configuration mean the indices of the cars and the empty cells, respectively. Thus 
the number of the cars K is 10 and the length of the circuit L is 38 in this case. We set the 
monitoring period no at 2. The speed of the cars 4, 9 and 0 is 3, which is the top speed no of 
this case. The speed of the car 5 is 2, whose headway is also 2. All the other cars’ speeds are 1, 
whose headways are also 1 except for the car 3. Thus the headways of the cars in tha past have 
nothing to do with the motion of the cars in the future except for the car 3. The headway of 
the car 3 at the time —1 is set to be 1. 

Out of the above initial configuration (13), the equation (1) generates flow of vehicles as 
follows 

0: I u 2u3uuu4uuuuuuu5uu6u7u8u9uuuuuuu0uuuu , 

1 ■ u1u2u3uuuuu4uuuuuu5u6u7lj8liulj9uuuuuuu0u , 

2: 0 u 1 u 2 u 3uuuuuuu4uuuu5u6u7u8uuuuu9uuuuuu, 

3: u0ulu2 u uu3uuuuuuu4uu5u6u7u8uuuuuuu9uuu, 

4: uu0ulu2uuuuu3uuuuuu4u5 u 6u7uuu8uuuuuuu9, 

5: u9u0u1u2|_iuuliljuu3ul]uu4u5lj6lj7|juuuu8uuuuu , 

6: uu9u0uluuu2uuuuuuu3uu4u5u6u7uuuuuuu8uu• 

Note that the minimum speed of the cars n“ in is 1 in the flow above. We notice that the 
configuration at the time 3 is obtained by moving all the cells of the initial configuration one 
cell rightward as well as changing the car indices k to k — 1 modulo 10. The configuration at the 
time 6 is also obtained by doing the same shifts and changes of car indices to the configuration 
at the time 3. In this sense, the above flow is a periodic motion of cars whose period is 3 in this 
case. The length of the jam, or the number of the cars running at the minimum speed, is thus 
almost constant. Intermediate speeds also appear but only temporarily. That is why we call 
them stationary flows of the s2s-OVCA. Roughly speaking, the slow flow we shall deal with is 
the stationary flow of the type shown above. The density of the cars p and the average flow Q 
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Figure 3. The slow flow of the s2s-OVCA. Cars in the cells drawn by dashed lines are omitted except 
for the car 1 at the time n + n o + 1. 


over the period, or the no + 1 = 3 steps, are calculated as 


P = 


K _ 5 
T ~ 19’ 


Q 


l 

(no + 1 )L 


K n 0 

EE»f 


k =0 n '=0 


6+3+3+5+9+4+3+3+3+9 
3 x 38 


8 

19’ 


which will be verified with the formula we shall derive shortly. 

Let us consider such slow flows as we have seen above as the specific solutions in a more 
general manner. Fig. 3 shows configurations of a slow flow at times n and n + no + 1. Since 
we employ the periodic boundary condition, two cells containing the car K are identified. As 
a property of the slow flow, we assume that the slow flow is periodic in the sense that the 
configuration at the time n + no + 1, which is shown in the box in Fig. 3, is given by the 
rightward displacement of the entire configuration at the time n in the box by nou“ in — 1 
cells. The flow provided by this displacement of the entire configuration in no + 1 time steps is 
—— p■ For example, the rightward displacement mentioned above for the slow flow in Fig. 3 
is 2x1 — 1 = 1, which agrees with the observation before. The set of stationary flows given 
in [15] has the property of the slow flow we here assume. 

Here we should note that the leftward displacement of the car K by L cells, namely whole 
the circuit length, which is fictitiously introduced to make the shifted initial configuration from 
the real configuration at the time no + 1 in the sense that the numerical order of the car arrays 
is maintained. In order to compensate the underestimation of the flow brought about by this 
leftward displacement, we have to add the flow corresponding to the rightward displacement of 
the car I\ by L cells in no + 1 time steps, („ (1 +i)l ’ ^ = no+ 1 ' Thus ^ ie fl° w of the slow flow 
with the minimum speed u+ n is given by 


n 0 + 1 


1 


n 0 + 1' 


0 < <° n < v 0 . 


(14) 


For example, substitution of p = no = 2 and w“ in = 1 into equation (14) yields 


_ 2x1—1 5 1 8 

Q = -x-1-= —, 

^ 2 + 1 19 2 + 1 19’ 


which agrees with the flow Q = for the slow flow given above as an specific solution. The 
formula (14) agrees with the flow-density relation given by numerical experiments, as we can 
see in Fig. 2. Three branches labeled with v = 2,1 and 0 are the flow-density relations with 
the minimum speeds v“ in = v in Fig. 2. Roughly speaking, the traffic that forms the branch 
corresponding to u“ in consists of groups of cars running at the top speed vo and other groups 
of cars running at the minimum speed u“ in , as the specific solution evolving from the initial 
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configuration (13) has shown. A set of specific solutions that correspond to the branches in the 
fundamental diagram is given in [15]. 

The maximum density faax( Cin ) that allows the minimum speed to be w“ in is 

Anax(^ m in) = T7' (15) 

The flow Q(/O m ax(t>mi n )) corresponding to the maximum density p m ax(^ mir< ) is then given by 

Q^max^min)) = Pmax(Unin)' u min- (15) 

Since the two equations (15) and (16) holds at the same time, they leads to Q(p m ax.(V m\ n )) + 
Pmax(^ min ) = 1- Thus all the end points of the branches must be on the line 

Q + P= 1. (17) 


The branching point, or the minimum density, of the flow-density relation of the slow flow 
corresponding to the minimum speed r* in is determined by the intersection of the flow density 
relations of the free flow and the slow flow 


Pmin 


oo \ 
min/ 


l 

M v o - Cl) + v 0 + 1' 


(18) 


In Fig. 2, the branching points corresponding to u“ in = 2,1 and 0 are encircled with small circles, 
which agree with the above formula (18). The density of the cars p needs to be sufficiently large 
so as to form the slow flow with the minimum speed 'u* in . The branching point gives the lower 
bound of such density. 

Since the s2s-OVCA (1) is a deterministic CA, the initial configuration determines the final 
state. Its numerical simulation is very robust against, or more precisely speaking, free from 
numerical errors. Thus all the stationary flows including the free flow and the slow flows 
beyond the minimum density Pmin(^mm) are stable in the numerical simulation. That is why 
we were able to obtain the flow-density relation with several “metastable” states consisted by 
the free flow and the slow flows beyond the minimum density, as was shown in Figs. 1 and 2. 
Though these metastable states are robust against numerical errors in simulation, but they are 
generally unstable against perturbation, which gives a reason of their name. For example, when 
one gives a perturbation to the headways that is equal to the minimum speed u“ in in the groups 
of cars running at the minimum speed u“ in , a car with a velocity that is less than appears. 
And such a perturbed car generally becomes a seed of a group of cars running at a velocity 
slower than i)^ in , which eventually slows down whole the traffic. That is why we call these 
states metastable, except for the slow line with u“ in = 0, which we cannot make slower. 

When the monitoring period no is zero, all of the slow flows (14) goes to the line of the end 
points of the metastable branches (17). Thus the monitoring period plays an essential role in 
the formation of the metastable states in the flow-density relation. Thus it could be determined 
by comparing the flow-density relation of the s2s-OVCA and observed ones. 


4 Summary 

We have shown an inverse ultradiscretization from the s2s-OVCA (1) to an integral-differential 
equation (9), which is an extension of the Newell model (10). Since the Newell model [8] and 
the s2s-OVCA [9] are extended models of the OV [1] and the s2s models [12] respectively, the 
s2s-OVCA is interpreted as a CA-type hybrid of the OV and the s2s models. 

Using the features of the stationary flows observed in the numerical experiments, we have 
derived the flow-density relations of the stationary flow of the s2s-OVCA. The flow-density 
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relations of the s2s-OVCA were numerically obtained [9] and then derived by use of a set of 
stationary flows [15]. 

Since the s2s-OVCA is a deterministic CA, the model is suitable for a simulation with a much 
bigger system size than the length of the circuit, L = 100, in our numerical simulation. We 
expected that it would be sufficient to capture the characteristics of stationary flows of the s2s- 
OVCA on the circuit, which is determined by the density of cars and initial configuration. In 
order to observe finite-size effects of open boundaries, for example, we expect that the s2s-OVCA 
will be a good tool. 

The s2s-OVCA has several types of monotonicity in its time evolution, which extend the 
results shown for the tlq = 1 case [13] . We expect that the nronotonicity determines the relaxation 
to the stationary flow from the initial configuration as well as the property of the stationary flow 
we assume here. We hope that results on the relaxation to stationary flows and the monotonicity 
in the time evolution of the s2s-OVCA will be reported soon. 
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